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Modeling the Spectral Energy Distributions and Variability of Blazars 
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In this review, recent progress in theoretical models for the broadband (radio through 7-ray) 
emission from blazars are summarized. The salient features of both leptonic and hadronic models 
are reviewed. I present sample modeling results of spectral energy distributions (SEDs) of different 
types of Fermi-detected blazars along the traditional blazar sequence, using both types of models. 
In many cases, the SEDs of high-frequency peaked blazars (HBLs) have been found to be well 
represented by simple synchrotron -I- synchrotron self-Compton (SSC) models. However, a few 
HBLs recently discovered as very-high-energy (VHE) gamma-ray emitters by VERITAS are actually 
better represented by either external- Compton or hadronic models. Often, spectral modeling with 
time-independent single-zone models alone is not sufficient to constrain models, as both leptonic 
and lepto-hadronic models are able to provide acceptable fits to the overall SED. This degeneracy 
can be lifted by considering further constraints from spectral variability. Recent developments of 
time-dependent and inhomogeneous blazar models will be discussed, including detailed numerical 
simulations as well as a semi-analytical approach to the time-dependent radiation signatures of 
shock-in-jet models. 



I. INTRODUCTION 

Blazars (BL Lac objects and 7-ray loud flat spec- 
trum radio quasars [FSRQs]) are the most extreme 
class of active galaxies known. Their broadband 
(from radio to 7-rays) continuum SEDs are dominated 
by non-thermal emission and consist of two distinct, 
broad components: a low-energy component from ra- 
dio through UV or X-rays, and a high-energy compo- 
nent from X-rays to 7-rays (see, e.g.. Figure [1]). 

Blazars are sub-divided into several types, de- 
fined by the location of the peak of the low-energy 
(synchrotron) SED component. Low-synchrotron- 
peaked (LSP) blazars, consisting of fiat-spectrum ra- 
dio quasars (FSRQs) and low-frequency peaked BL 
Lac objects (LBLs), have their synchrotron peak in 
the infrared regime, at Vg ^ 10^^ Hz. Intermediate- 
Synchrotron-Peaked (ISP) blazars, including LBLs 
and intermediate BL Lac objects (IBLs), have their 
synchrotron peak at optical - near-UV frequencies at 
lO^'^Hz < Vs < 10^^ Hz, while High-Synchrotron- 
Peaked (HSP) blazars, which are almost all high- 
frequency-peaked BL Lac objects (HBL), have their 
synchrotron peak at X-ray energies with Vs > Hz 
Q . Such a sequence appears to be more physical than 
the historical distinction between quasars and BL Lac 
objects based on the equivalent width of optical emis- 
sion lines [5^, as the equivalent width is dependent 
on the strongly time- variable flux of the non-thermal 
continuum [5l|. 

The blazar sequence had first been identified by [1^ , 
and associated also with a trend of overall decreas- 
ing bolometric luminosity as well as decreasing 7-ray 
dominance along the sequence FSRQ LBL HBL. 
According to this classification, the bolometric power 
output of FSRQs is strongly 7-ray dominated, espe- 
cially during flaring states, while HBLs are expected 
to be always synchrotron dominated. However, while 



the overall bolometric-luminosity trend still seems to 
hold, recently, even HBLs seem to undergo episodes 
of strong 7-ray dominance. A prominent example was 
the giant VHE 7-ray flare of the HBL PKS 2155-304 
in 2006, during which the SED was strongly 7-ray 
dominated [12]. 

Figure [1] shows examples of blazar SEDs along the 
blazar sequence, from the FSRQ 3C279 (a), via the 
LBL BL Lacertae (b) and the IBL 3C 66A (c), to the 
HBL RGB J0710-f 591 (d). The sequence of increasing 
synchrotron peak frequency is clearly visible. How- 
ever, the Fermi spectrum of the LBL BL Lacertae 
indicates a 7-ray flux clearly below the synchrotron 
level, while the SED of the IBL 3C 66A is clearly 
dominated by the Fermi 7-ray flux, in contradiction 
with the traditional blazar sequence. 

Blazars are known to be variable at all wavelengths. 
Typically, the variability amplitudes are the largest 
and variability time scales are the shortest at the 
high-frequency ends of the two SED components. In 
HBLs, this refers to the X-ray and VHE 7-ray regimes. 
High-energy variability in some HBLs (PKS 2155- 
305, Mrk 501, Mrk 421) has been observed on time 
scales down to just a few minutes [l^, [13, SII. The 
flux variability of blazars is often accompanied by 
spectral changes which are sometimes associated with 
inter-band or intra-band time lags. In a few HBLs 
(especially, Mrk 421), variability patterns have been 
found which can be characterized as spectral hystere- 
sis in hardness-intensity diagrams [e.g. .1431 [571 [85l[91|. 
However, even within the same object this feature 
tends not to be persistent over multiple observations. 
Also in other types of blazars, hints of time lags be- 
tween different observing bands are occasionally found 
in individual observing campaigns (e.g., [l^Hj]), but 
the search for time-lag patterns persisting through- 
out multiple years has so far remained unsuccessful 
[see, e.g.,[5j, for a systematic search for time lags be- 
tween optical, X-ray and 7-ray emission in the quasar 
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FIG. 1: Spectral energy distributions of four sub-classes of blazars: a) the FSRQ 3C279 (from ;3Q]), b) the LBL BL 
Lacertae, (data from Q]), c) the intermediate BL Lac 3C66A (data from Q]), and d) the HBL RGB J0710-I-591 (data 
from [^). In Panel a) (3C279), lines are one-zone leptonic model fits to SEDs at various epochs shown in the figure. 
In all other panels, red lines are fits with a leptonic one-zone model; green lines are fits with a one-zone lepto-hadronic 
model. 



3C279]. 



II. LEPTONIC AND HADRONIC BLAZAR 
MODELS 

The high apparent bolometric luminosities (assum- 
ing isotropic emission) , rapid variability, and apparent 
superluminal motions of radio components in blazar 
jets provide compelling evidence that the nonthermal 
continuum emission of blazars is produced in <, 1 light 
day sized emission regions, propagating relativistically 
with velocity /3rc along a jet directed at a small angle 
^obs with respect to our line of sight [e.g., [t^. Let 
r = (1 — /3r)~^/^ be the bulk Lorentz factor of the 
emission region, then Dopplcr boosting is determined 
by the Doppler factor D = (r[l — /3r cos^obs])^^- 
In the following, primes denote quantities in the co- 
moving frame of the emission region. The observed 



frequency v^hs is related to the emitted frequency 
through Vo\is = Dv' /{I -\- z), where z is the red- 
shift of the source. The energy fluxes are connected 
through F°^^^ = F^, . Intrinsic variability on a co- 
moving time scale t'^^^ will be observed on a time scale 
^var = ^var (1 + z) / D . Using the latter transformation 
along with causality arguments, any observed vari- 
ability leads to an upper limit on the size scale of the 
emission region through R<ct°^^ + z). 

It is generally agreed that the low-frequency compo- 
nent in the blazar SED arises from synchrotron emis- 
sion of relativistic electrons. However, there are two 
fundamentally different approaches concerning the 
high-energy emission. If protons are not accelerated 
to ultrarelativistic energies (exceeding the threshold 
for p7 pion production on synchrotron and/or exter- 
nal photons), the high-energy radiation will be domi- 
nated by Compton emission from ultrarelativistic elec- 
trons and/or pairs (leptonic models). In the oppo- 
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site case, the high-energy emission will be dominated 
by cascades initiated by p7 pair and pion produc- 
tion as well as proton, tt^, and /i^ synchrotron ra- 
diation, while primary leptons are still responsible 
for the low-frequency synchrotron emission (hadronic 
models). The following sub-sections provide a brief 
overview of the main radiation physics aspects of both 
leptonic and hadronic models. A more in-depth re- 
view of the radiation physics in relativistic jets can be 
found in (ssj. 

When fitting SEDs of blazars, consideration should 
be given to the energy balance between the magnetic 
field and the particle content in the jet, as this con- 
tains information on the jet launching and accelera- 
tion mechanisms. If the relativistic jets of AGN are 
powered by the rotational energy of the central black 
hole [1^ , the jets are expected to be initially Poynting- 
flux dominated. The energy thus carried primarily in 
the form of magnetic fields then needs to transferred 
to relativistic particles which can then produce the 
observed high-energy emission. This energy conver- 
sion is expected to cease when the energy densities 
in the magnetic field and in relativistic particles ap- 
proach equipartition. In this scenario, the jets are 
therefore not expected to become matter dominated 
within the central few parsecs of AGN, where the high- 
energy emission in blazars is believed to be produced 
[ssj . Also, if magnetic pressure plays an essential role 
in collimating and confining the jets out to pc or kpc 
scales, the energy density (and hence pressure) of par- 
ticles in the jet can not dominate over the magnetic- 
field pressure since otherwise the jets would simply 
expand conically. 

If SED fits require far sub-equipartition magnetic 
fields, this may indicate that the jets are magneto- 
hydrodynamically powered by the accretion flow. In 
this case, magnetic fields may be self-generated in 
shear fiows in the case of a fast inner spine sur- 
rounded by a slower, mildly relativistic outer fiow 
(sheath) . Kelvin-Helmholtz type instabilities can then 
lead to anisotropic particle acceleration, and the self- 
generated magnetic fields are expected to remain be- 
low equi part ition with the thermal energy carried in 
the flow [ll[6l|. 

A. Leptonic models 

In leptonic models, the high-energy emission is pro- 
duced via Compton upscattering of soft photons off 
the same ultrarelativistic electrons which are produc- 
ing the synchrotron emission. Both the synchrotron 
photons produced within the jet [SSC = Synchrotron 
Self-Compton: l25l l66ll67|. and external photons (EC 
= External Compton) can serve as target photons 
for Compton scattering. Possible sources of exter- 
nal seed photons include the accretion disk radiation 
[e.g., ,37^ ,38||, reprocessed optical - UV emission from 



circumnuclear material [e.g., the BLR: [22l. [sol l46l. [soj , 

infrared emission from warm dust , or synchrotron 
emission from other (faster/slower) regions of the jet 
itself gl nil. 

If the observed blazar emission were arising in a 
stationary emission region (whose size is constrained 
by causality arguments), its compactness would, in 
many cases, be too large for high-energy 7-rays to es- 
cape from it. This problem is overcome by relativistic 
Doppler boosting so that model fits can be achieved 
with parameters such that the 77 opacity of the emis- 
sion region is low throughout most of the high-energy 
spectrum. However, at the highest photon energies, 
this effect may make a non-negligible contribution to 
the formation of the emerging spectrum [TT| and lead 
to the re-processing of some of the radiated power 
towards lower frequencies. 77 absorption within the 
AGN may also occur on photons external to the emis- 
sion region, as required to be present if 7-rays are 
produced through the EC mechanism [e.g.l62l.l74.[76|. 
This absorption process might lead to Compton sup- 
ported pair cascades which are expected to release 
their peak power at MeV - GeV energies [e.g..[77l,[78|. 

Also the deceleration of the jets may have a sig- 
nificant impact on the observable properties of blazar 
emission through the radiative interaction of emission 
regions with different speed [3, . The deceleration 
of the high-energy emission region is expected to cause 
characteristic variability features which may be remi- 
niscent of those seen in gamma-ray burst light curves 
[32} . Varying Doppler factors may also be a result of a 
slight change in the jet orientation without a substan- 
tial change in speed, e.g., in a helical-jet configuration 
[e.g., nil- In the case of ordered magnetic- field struc- 
tures in the emission region, such a helical configura- 
tion should have observable synchrotron polarization 
signatures, such as the prominent polarization-angle 
swing recently observed in conjunction with an opti- 
cal + Fermi 7-ray flare of 3C 279 Q. 

In the simplest models, the underlying lepton (elec- 
trons and/or positrons) distribution is ad-hoc pre- 
specified, either as a single or broken power-law with 
a low- and high-energy cut-off. While leptonic mod- 
els under this assumption have been successful in 
modeling blazar SEDs [e.g., |47|, they lack a self- 
consistent basis for the shape of the electron distri- 
bution. A more realistic approach consists of the self- 
consistent steady-state solution of the Fokker-Planck 
equation governing the evolution of the electron dis- 
tribution function, including a physically motivated 
(e.g., from shock-acceleration theory) acceleration of 
particles and all relevant radiative and adiabatic cool- 
ing mechanisms [e.g., H, HO, 113 ■ The model described 
in S, b ased on the time-dependent model of [2^ [see 
also [34| . has been used to produce the leptonic model 
fits shown in Fig. [TJ 

In order to reproduce not only broadband SEDs, 
but also variability patterns, the time-dependent elec- 
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tron dynamics and radiation transfer problem has to 
be solved self-consistently. Such time-dependent SSC 
models have been developed by, e.g., [13, 113, Issl. l82j|. 
External radiation fields have been included in such 
treatments in, e.g., [26l. IsTI. Issj . 

Leptonic models have generally been very success- 
fully applied to model the SEDs and (in few cases) 
spectral variability of blazars. The radiative cooling 
time scales (in the observers's frame) of synchrotron- 
emitting electrons in a typical J5 ~ 1 G magnetic field 
are of the order of several hours - ~ 1 d at optical 
frequencies and ;<! hr in X-rays and hence compati- 
ble with the observed intra-day variability. However, 
the recent observation of extremely rapid VHE 7-ray 
variability on time scales of a few minutes poses severe 
problems to simple one-zone leptonic emission models. 
Even with large bulk Lorentz factors of ^ 50, causal- 
ity requires a size of the emitting region that might 
be smaller than the Schwarzschild radius of the central 
black hole of the AGN [2l|. As a possible solution, it 



has been suggested [87| that the 7-ray emission region 
may, in fact, be only a small spine of ultrarelativistic 
plasma within a larger, slower-moving jet. Such fast- 
moving small-scale jets could plausibly be powered by 
magnetic reconnection in a Poynting-flux dominated 
jet, as proposed by [4^. 



B. Hadronic models 

If a significant fraction of the jet power is con- 
verted into the acceleration of relativistic protons 
in a strongly magnetized environment, reaching 
the threshold for pj pion production, synchrotron- 
supported pair cascades will develop, initiated by pri- 
mary 7r° decay photons and synchrotron emission from 
secondary pions, muons and electrons/positrons at 
ultra-high 7-ray energies, where the emission region 
is highly opaque to 77 absorption [H, . The ac- 
celeration of protons to the necessary ultrarelativistic 
energies (i?™'*'^> 10^^ eV) requires high magnetic fields 
of several tens of Gauss to constrain the Larmor radius 
i?L = 3.3 X lO^'^ Bf^ Ei9 cm, where B = 10 Bi G, and 
Ep — 10^^ Eig eV, to be smaller than the size of the 
emission region (typically inferred to be i?<,10^^ cm 
from the observed variability time scale) . In the pres- 
ence of such high magnetic fields, the synchrotron 

111 and of sec- 
must be taken 




radiation of the primary protons 
ondary muons and mesons (tiI - It^ 
into account in order to construct a self-consistent 
synchrotron-proton blazar (SPB) model. Electromag- 
netic cascades can be initiated by photons from tt^- 
decay ("tt" cascade"), electrons from the tt^ — ^ — > 
decay ( "tt^ cascade" ) , p-synchrotron photons ( "p- 
synchrotron cascade"), and /Lt-, tt- and X-synchrotron 
photons ( "/z^-synchrotron cascade"). 

(7^ and [t^] have shown that the "tt^ cascades" 
and "tt^ cascades" from ultra-high energy protons 



generate featureless 7-ray spectra, in contrast to "p- 
synchrotron cascades" and "/x^ -synchrotron cascades" 
that produce a two-component 7-ray spectrum. In 
general, direct proton and fi^ synchrotron radiation 
is mainly responsible for the high energy bump in 
blazars, whereas the low energy bump is dominanted 
by synchrotron radiation from the primary e~, with a 
contribution from secondary electrons. 

Hadronic blazar models have so far been very dif- 
ficult to investigate in a time-dependent way be- 
cause of the very time-consuming nature of the re- 
quired Monte-Carlo cascade simulations. In general, 
it appears that it is difficult to reconcile very rapid 
high-energy variability with the radiative cooling time 
scales of protons, e.g., due to synchrotron emission, 
which is t°^^ = 4.5 X 10^ (1 + z)D:[^B^'^ E^^ s [§], 
i.e., of the order of several days for ^ 10 G magnetic 
fields and typical Doppler factors D = lOZ^i. How- 
ever, rapid variability on time scales shorter than the 
proton cooling time scale may be caused by geomet- 
rical effects, for example in a helical jet. 

In order to avoid time-consuming Monte-Carlo sim- 
ulations of the hadronic processes and cascades in- 
volved in lepto-hadronic models, one may consider a 
simplified prescription of the hadronic processes, (ssj l 
have produced analytic fit functions to Monte-Carlo 
generated results of hadronic interactions using the 
SOFIA code [t^I. Those fits describe the spectra of 
the final decay products, such as electrons, positrons, 
neutrinos, and photons from tt*^ decay. This neglects 
the possibility of substantial synchrotron losses to 
muons and pions before they decay. Since the muon 
life time is longer than the pion life time, the more 
stringent condition is given by the requirement that 
isy.^i. S> td,^n where igy and denote the synchrotron 
and decay time scales, respectively. This can be trans- 
lated into the requirement that B <^ 5-67p 10 
where 7p,max,io is the Lorentz factor of the highest- 
energy protons in units of 10^°. [H^ have developed 
analytical fits to the secondary-production spectra in 
hadronic processes for all individual steps, thus al- 
lowing one to take muon and pion synchrotron (and 
Compton) emission into account properly. 

A hadronic blazar model based on the templates of 
[58| has been developed by [3J] [see alsoM. This 
model employs a semi-analytical solution for UHE 
7-ray induced, synchrotron-supported pair cascades. 
The hadronic fits shown in Figfl] b) - d) have been 
performed with this model. As the low-frequency 
component is electron-synchrotron emission from pri- 
mary electrons, it is not surprising that virtually iden- 
tical fits to the synchrotron component can be pro- 
vided in both types of models. In the high-frequency 
component, strongly peaked spectral shapes, as, e.g., 
in 3C 66A and RGB J0710+591 require a strong 
proton-synchrotron dominance with the cascading of 
higher-energy (3> TeV) emission only making a mi- 
nor contribution to the high-energy emission. This, 
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in fact, makes it difficult to achieve a substantial ex- 
tension of the escaping high-energy emission into the 
> 100 GeV VHE 7-ray regime. In objects with a 
smoother high-energy SED, e.g., BL Lacertae in Fig. 
[TJj, a substantially larger contribution from cascade 
emission (and leptonic SSC emission) is allowed to 
account for a relatively high level of hard X-ray / soft 
7-ray emission. This also allows for a substantial ex- 
tension of the 7-ray spectrum into the VHE regime. 



III. SED FITS ALONG THE BLAZAR 
SEQUENCE 

In the framework of leptonic models, the blazar se- 
quence LSP ISP — > HSP is often modeled through 
a decreasing contribution of external radiation fields 
to radiative cooling of electrons and associated EC 
7-ray emission ^7|. In this sense, HBLs have tradi- 
tionally been well represented by pure synchrotron- 
self-Compton models, while FSRQs usually require 
a substantial EC component. This interpretation is 
consistent with the observed strong emission lines in 
FSRQs, which are absent in BL Lac objects. At the 
same time, the denser circumnuclear environment in 
quasars might also lead to a higher accretion rate and 
hence a more powerful jet, consistent with the over- 
all trend of bolometric luminosities along the blazar 
sequence. This may even be related to an evolution- 
ary sequence from FSRQs to HBLs governed by the 
gra dual depletion of the circumnuclear environment 

However, in this interpretation, it would be ex- 
pected that mostly HBLs (and maybe IBLs) should 
be detectable as emitters of VHE 7-rays since in LBLs 
and FSRQs, electrons are not expected to reach ~ TeV 
energies. This paradigm has been undermined by 
the recent VHE 7-ray detections of lower-frequency 
peaked objects such as W Comae [1], 3C66A 
PKS 1424-H240 0, BL Lacertae fTl, 85 0716+714 
2a, and even the FSRQs 3C 279 jll|, PKS 1222+21 
H and PKS 1510-089 

The overall SEDs of IBLs detected by VERITAS 
could still be fit satisfactorily with a purely leptonic 
model. Fitting the SEDs of the IBLs 3C66A and 
W Comae with a pure SSC model, while formally 
possible, would require rather extreme parameters. 
In particular, magnetic fields several orders of mag- 
nitude below equipartition would be needed, which 
might pose a severe problem for jet collimation (see 
the discussion in Section |n|. Much more natural fit 
parameters can be adopted when including an EC 
component with an infrared radiation field as target 
photons d, @. 

In j2^ it has been demonstrated that the VHE 7-ray 
detection of the FSRQ 3C 279 poses severe problems 
for any variation of a single-zone leptonic jet model. 
The fundamental problem for the leptonic model in- 



terpretation is the extremely wide separation between 
the synchrotron and 7-ray peak frequencies. In a 
single-zone SSC interpretation this would require a 
very high Lorentz factor of electrons at the peak of 
the electron distribution and, in turn, an extremely 
low magnetic field. In an EC interpretation the SSC 
component, usually dominating the X-ray emission in 
leptonic fits to FSRQs like 3C 279, is too far sup- 
pressed to model the simultaneous RXTE spectrum. 
As an alternative, the optical emission could be pro- 
duced in a separate emission component. A pure SSC 
fit to the X-ray and 7-ray component is technically 
possible, but also requires a far sub-equipartition mag- 
netic field. Much more natural parameters could be 
achieved in a fit with the synchrotron-proton-blazar 
model of 

A. HBLS disfavouring SSC 

Until very recently, the SEDs of essentially all HBLs 
detected in VHE 7-rays could be well modelled us- 
ing simple leptonic SSC models. The most attrac- 
tive feature of this model is its relative simplicity, 
requiring only a small number of fit parameters. If 
the peak frequencies and corresponding fluxes in a 
blazar SED are well measured and an additional con- 
straint on the emission region size can be obtained 
from rapid variability, then all relevant parameters 
can be uniquely determined [e.g.,[8^. However, recent 
new discoveries by VERITAS with good simultaneous 
or quasi-simultaneous multiwavelength coverage pro- 
vided strongly constraining SEDs that were difficult 
to reproduce with such a model. 

The HBL RX J0648. 7+1516 (z = 0.179) was de- 
tected as a VHE7-ray source by VERITAS in 2010 
March - April [131 . Simultaneous multiwavelength 
coverage was provided by Fermi-LAT, Swift XRT and 
UVOT and the MDM Observatory (optical B, V, R 
bands). The resulting SED is shown in the left panel 
of Figure O An attempt to fit this SED with a 
pure SSC model clearly overshoots both the simul- 
taneous Fermi upper limits and the 2-year average 
spectrum. This is a consequence of the apparent sharp 
spectral break between the Fermi and VERITAS en- 
ergy ranges, which is not consistent with the smooth 
curvature of an SSC spectrum (resulting from the 
broad ranges of scattering electron energies and syn- 
chrotron seed photon frequencies, in combination with 
the gradual onset of Klein- Nishina effects). A much 
more satisfactory leptonic fit can be achieved when al- 
lowing for an EC contribution to the 7-ray emission. 
In the fit shown by the dot-dashed curve in the left 
panel of Fig. [51 an external blackbody radiation field 
at a temperature of T = 1000 K has been assumed, 
which could result from an infrared-emitting dust re- 
gion around the AGN. Its luminosity corresponds to 
a fiux of vF^ <l{f Jy Hz, in agreement with its non- 
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FIG. 2: Spectral energy distributions of the HBLs RX J0648. 7+1516 [left;!?! and RBS 0413 [right;!!!]. The Fermi upper 
limits for RX J0648. 7+1516 correspond to the ~ 2 months time window simultaneous with the VERITAS observations, 
while the brown Fermi points show the 2-year average spectrum. The solid curve indicates a pure leptonic SSC fit. The 
dashed dot-dashed curve is a leptonic SSC + EC fit, while the dashed curve shows a fit with the semi-analytical hadronic 
model described in Section Til Bl See text for details of the individual fits. The blue curve in the right panel represents a 
phenomenological representation of the (thermal dominated) host-galajcy emission. 



detection in the SED. 

Alternatively, the hadronic model described above 
also provides a satisfactory fit to the SED (the dashed 
curve in the left panel of Fig. [2]). The model requires 
protons to be accelerated up to Sp^max = l.SxlO^^eV, 
with a total kinetic luminosity in relativistic protons 
of Lp — 2 X 10'*^ erg/s. While hadronic models often 
suffer from extreme power requirements, this appears 
to be a very reasonable jet power. The magnetic field 
of i? = 10 G is slightly above equipartition with the 
energy content in relativistic protons. 

The HBL RBS 0413 {z = 0.19) was detected by 
VERITAS in data accumulated between 2008 Septem- 
ber and 2010 January. Contemporaneous multiwave- 
length coverage was provided by Fermi-LAT, Swift 
XRT and UVOT, and MDM (R-band). A leptonic 
SSC model fit (solid curve in the right panel of Fig. 
[2]) produces a GeV 7-ray spectrum significantly harder 
than indicated by Fermi-LAT, and requires a mag- 
netic field energy density a factor ~ 6 x 10~^ below 
equipartition. A much more satisfactory represen- 
tation of the entire SED (including the Fermi-LAT 
spectrum), can be achieved by including an EC com- 
ponent with a thermal blackbody radiation field at 
T = 1.5 X 10^ K. This is shown by the dot-dashed 
curve in the right panel of Figure . Parameters in 
almost exact equipartition between the magnetic field 
and the relativistic electron population can be chosen, 
and the postulated infrared radiation field corresponds 
to a fiux of i/Fl^ ~ 5 X 10® JyHz, consistent with the 
absence of direct evidence for it in the infrared spec- 
trum of the source. 

The dashed curve in the right panel of Fig. [5] shows 
a hadronic fit to the SED of RBS 0413, which also 



provides a better representation than the pure SSC 
model. For this fit, protons need to be accelerated 
to iJp^max = 1.6 X 10^^ eV with a kinetic power of 
= 2 X 10''*^ erg/s. The magnetic field of B = 30 G is 
in almost exact equipartition with the energy content 
in relativistic protons. 

Thus, both RX J0648. 7+1516 and RBS 0413 ex- 
hibit SEDs that disfavor a pure leptonic single-zone 
SSC interpretation, while both an EC model using a 
thermal infrared radiation field, and a hadronic model 
provide convincing fits. Both EC and hadronic model 
fits can be achieved with with reasonable energy re- 
quirements and parameters close to equipartition be- 
tween the magnetic field and the dominant relativistic 
particle population. 



IV. INHOMOGENEOUS JET MODELS 

The complicated and often inconsistent variability 
features found in blazars provide a strong motivation 
to investigate jet models beyond a simple, spherical, 
one-zone geometry. The idea behind phenomenologi- 
cal multi-zone models like the spine-sheath model of 
[S?! or the decelerating-jet model of Q was that dif- 
ferential relativistic motion between various emission 
zones will lead to Doppler boosting of one zone's emis- 
sion into the rest frame of another zone. This can re- 
duce the requirements of extreme bulk Lorentz (and 
Doppler) factors inferred from simple one-zone lep- 
tonic modeling of rapidly variable VHE 7-ray blazars, 
such as Mrk 501 or PKS 2155-304, and has led to suc- 
cessful model fits to the SEDs of those sources with 
much more reasonable model parameters. 
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FIG. 3: Dependence of the SED characteristics of time- 
averaged spectra from the intermal-shock model, on the 
external radiation energy density iicxt- 



However, the models mentioned above do not treat 
the radiation transport and electron dynamics in a 
time-dependent way and do therefore not make any 
robust predictions concerning variability and inter- 
band cross-correlations and time lags. In order to 
address those issues, much work has recently been de- 
voted to the investigation of the radiative and timing 
signatures of shock-in-jet models, which will be sum- 
marized in the following sub-section. 



A. Shock-in-jet models 

Early versions of shock-in-jet models were devel- 
oped with focus on explaining radio spectra of ex- 
tragalactic jets, e.g., by Their application to 

high-energy spectra of blazars was proposed by [83 |. 
Detailed treatments of the electron energization and 
dynamics and the radiation transfer in a standing 
shock (Mach disk) in a blazar jet were developed by 
[13,182, 83] . The internal-shock model discussed in [69| 
and 56] assumes that the central engine is intermit- 
tently ejecting shells of relativistic plasma at varying 
speeds, which subsequently collide. Such models have 
had remarkable success in explaining SEDs and time 
lag features of generic blazars. 

The realistic treatment of radiation transfer in an 
internal-shock model for a blazar requires the time- 
dependent evaluation of retarded radiation fields orig- 
inating from all parts of the shocked regions of the 
jet. The model system is therefore highly non-linear 
and can generally only be solved using numerical sim- 
ulations (e.g., [53, [63, Is^l)- As the current detailed 
internal-shock models employ either full expressions 
or accurate approximations to the full emissivities of 
synchrotron and Compton emission, a complete simu- 
lation of the time-dependent spectra and light curves 



is time-consuming and does therefore generally not 
allow to efhciently explore a large parameter space. 
General patterns of the SED, light curves and ex- 
pected time lags between different wavelength bands 
have been demonstrated for very specific, but obser- 
vationally very poorly constrained, sets of parame- 
ters. The most detailed treatment of the radiation- 
transfer aspect in the context of this models, includ- 
ing all light-travel time effects, has been presented in 
[35| , who employ a time-dependent multi-zone Monte- 
Carlo (for the radiation transfer) and Fokker-Planck 
(for electron dynamics) code. 

In order to be able to explore timing signatures 
of the internal-shock model for blazars on the basis 
of a broad parameter study, [3l| have developed a 
semi-analytical scheme for this model. The time- and 
space-dependent electron spectra, affected by shock 
acceleration behind the forward and reverse shocks, 
and subsequent radiative cooling, are calculated fully 
analytically. Taking into account all light travel 
time effects, the observed synchrotron and external- 
Compton spectra are also evaluated fully analytically, 
using a (5-function approximation to the emissivities. 
The evaluation of the SSC emission still requires a 
two-dimensional numerical integration. 

This semi- analytical model allowed the authors to 
efficiently scan a substantial region of parameter space 
and discuss the dependence on the characteristics of 
time-averaged SEDs, as well as cross-band correla- 
tions and time lags. As an example, Fig. [3] shows the 
dependence of the SED characteristics on the external 
radiation field energy density, Uoxt- In the classical 
interpretation of the blazar sequence, an increasing 
Ucxt corresponds to a transition from BL Lac spec- 
tral characteristics to FSRQ-like characteristics. The 
decreasing synchrotron peak frequency and increasing 
Compton dominance found in the parameter study are 
reproducing this effect. 

From the internal-shock simulations, energy- 
dependent light curves could be extracted. Usin g th e 
standard Discrete Correlation Function (DCF: |40| ) 
analysis, inter-band time lags could be extracted for 
any set of parameters. Fig. |4]shows the dependence of 
the inter-band time lags between optical. X-ray and 
Fermi 7-ray light curves (top panel) and the qual- 
ity of the cross correlation, characterized by the peak 
value of the DCF (bottom panel) as a function of the 
external radiation energy density. 

One of the most remarkable results of this study 
was that only slight changes in physical parameters 
can lead to substantial changes of the inter-band time 
lags and even a reversal of the sign of the lags. This 
may explain the lack of consistency of lags even within 
the same source. 
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